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' In order to study non-zero spin excitations of the recently proposed a- cluster conden- 

sation in the self-conjugate 4n nuclei, spatial deformation is introduced into the model wave 

. function of the a-cluster condensate. The rotational band states of *Be are investigated as 

^ ' a first step of a test case for the study of the deformation of the a-cluster condensate. Cal- 

fS) I culations reproduce well the binding energy of the 0"'' ground state and also the excitation 

^SJ ■ energy of the 2"*" state. Our 0"*" wave function is found to be almost exactly equal to the 0^ 

■ wave function obtained by the generator coordinate method using Brink's 2a wave function, 
f — ' The study shows that both the O'*' ground and 2^ excited states can be considered as having 

I a gas-like (i.e. weakly bound) 2a-cluster structure. As for the 0''' ground state, the change 

fSJ . of the wave function, due to the introduction of the deformation, is found to be very small. 

O ; 

I §1. Introduction 

^ . Recently we presented a conjecture that near the na threshold in self-conjugate 

pH I 4n nuclei there exist excited states of dilute density which are composed of a weakly 

^ ' interacting gas of alpha particles and which can be considered as an na condensed 

. stated. This conjecture was examined in ^^C and ^^O by using a new a-cluster wave 

r> I function of the a-particle condensate type. The second 0"*" state in ^^C near the 3a 

• threshold was interpreted as such a condensed state as well as the fifth 0^ T=0 state 

in ^^O near the 4a-particle threshold. 

The new a-cluster wave function of the a-particle condensate type which we used 
in our previous study is an eigen state of zero angular momentum and represents 
a condensation of a-particles in a spherically symmetric state. Therefore this new 
a-cluster wave function can describe only states of zero angular momentum. On 
the other hand, in ^-^C there exist some excited states above the 3a threshold which 
have non-zero angular momenta and are considered to have well-developed 3a cluster 
structure Actually many microscopic 3a cluster model calculations gave the result 
that not only the second 0"*" state near the 3a threshold but also some excited states 
with non-zero spin above the 3a threshold have 3a structure of dilute density u\ The 
question then arises whether we can identify some excitation modes of the a-cluster 
condensate with those excited states. 

As a possible excitation mode of the a-cluster condensate we can of course 
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consider partial breaking of the condensate like that due to the promotion of one 
a-cluster from the condensed s-orbit to a higher-lying orbit with non-zero spin. As 
another possibility for the excitation mode we can consider spatial deformation of the 
a-cluster condensate. If the deformation is not energetically stable, we can expect 
vibrational excitations around the spherical shape, while if the deformation is stable, 
we can expect the appearance of rotational excitations. 

The purpose of the present paper is to study the possibility of spatial deformation 
of the a-cluster condensate. This study is made in introducing a wave function of the 
a-cluster condensate with deformation which is obtained by a natural modification of 
our previous wave function of the spherical a-cluster condensate. As a first step for 
the study of deformation, we investigate the ^Be nucleus. The ground state of ^Be 
is known to have a well-developed 2a cluster structure with large spatial separation 
between the two a's in relative s-state. With only two a-particles this state can 
hardly be considered as an a-particle condensation, but rather as a weakly bound 
state of two bosons in relative s-state thus exhibiting a gas-like structure of dilute 
density. Our wave function of Ref.cP is perfectly adapted to also treat only two a's 
and it will serve us as a test case to treat more a-particles in the future. The nucleus 
®Be is known to have rotational excited states with large a-decay widths Fa, namely 
a 2+ state at 2.90 MeV with F^ = 1.45 MeV and a 4+ state at 11.4 MeV with ~ 7 
MeVo'. We will see that our investigation shows that in ^Be the a-cluster condensed 
state is deformable and the experimentally observed rotational band can be explained 
as being generated from the a-cluster condensate with deformation. We will also see 
that our 0"*" wave function is almost exactly equal to the 0"*" wave function obtained 
by the generator coordinate method using Brink's 2a wave function. 

§2. Deformed a-particle condensate 

2.1. Wave function 

The wave function of the a-cluster condensate with deformation can be expressed 
as follows by a slight modification of our previous wave function of the spherical a- 
cluster condensation: 

l^na) = (Ci)"|vac), (2-1) 

where the a-particle creation operator in a deformed center of mass orbit is given 
by 



Here 




and oJ.^(r) is the creation operator of a nucleon with spin-isospin err at the spatial 
point r. Our previous wave function of the spherical a-cluster condensation is just 
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the case of Px = Py = Pz = Ro in Eq.(2.2). The total na wave function therefore 
can be written as 

• r 4^1(7 4nT4n\^ no) 

'2X1 2X2, 2X 



oc A 



{ - + ^ + ^) } -^("l) • • • 4>{an)] , (2-4) 



where B'^ = + 2/3|, {k = x,y,z), and Xi = (l/4)X]n'"m is the center-of-mass 
coordinate of the ith a-cluster a,. The internal wave function of the a-cluster is 
(j){ai) oc exp[— (1/262) J2ni''^in — ^i)^]- The operator A is the total antisymmetrizer. 
It is to be noted that the wave function of Eqs. (2.1) and (2.4) expresses the state 
where n a-clusters occupy the same deformed center of mass orbit cxp{ — (-^X^ + 

-^Xy + -^Xl)^, while the intrinsic a-particle wave function stays spherical. We can 

easily see that the total center-of-mass motion is separated out of the wave function 
of Eqs. (2.1) and (2.4) because of the relation, 

n n n 

J2XI = nX^Gk + T^i^ik - Xokf, {k = X, y, z), Xg = (1/n) ^ X,. (2-5) 

1=1 i=l i=l 

The wave function of the a-cluster condensation with deformation expressed 
by Eqs. (2.1) and (2.4) is a superposition of Brink's a-cluster model wave function, 
^i^(i?i,---,i?„), 



na 



^na{Px,Py,Pz) = {riam, ■ ■ ■ r4nCr4nUn\^ 

= fd^Ri---d'Rn exp{-Y,[^ + ^ + ^)}^^{Rr,---,Rn), (2-6) 

i—i Px Py Pz 
^^{Rl, ■■■ , Rn) = det{(posiri - Rl)Xa-,Ti ■ ■ ■ fOs{r4n - Rn)Xai„Ti„}, (2-7) 

where Xar is the nucleon spin-isospin wave function. The form of Eq.(2.6) is directly 
obtained from Eqs. (2.1) ~ (2.3), and the form of Eq.(2.4) can be obtained from 
Eq.(2.6) using the relation, 

:det{(/5os(ri - R)x 



Vi! 

= (-l)'/%xp{ - ^{X - Rf}cj^{a), (2-8) 

where X = (1/4) ^^^^rj. An a-condcnscd wave function with good angular mo- 
mentum is obtained by projecting out the angular momentum from the deformed 
a-condensed wave function as 

^i^{P„Py,P,)= J df2Di^K{^)R{f2)^na{Px,Py,Pz) 

= j dQDi^K{Q) j d^Ri ■■■d^Rn exp(-^ ^ ^) 

i=l k=x,y,z 
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= J df2Di*^{f2) J d'R, ...d^i^ exp{ E ^^^^^^} 

i=l k=x,y,z k 

x^^{Ri,---,Rn), (2-9) 

where i? is the Euler angle, D'l^j^{f2) is the Wigner D-function, R{Q) is the rotation 
operator, and R{Q) is the 3x3 rotation matrix corresponding to R{Q). Here use 
is made of the relation, ipQs{R{^^)r — R) = ipos{f — R''^{f2)R). It is to be noted 
that the rotation with respect to the total angular momentum is equivalent to the 
rotation with respect to the total orbital angular momentum, since the intrinsic spins 
of a-clusters are saturated. 

In this paper we only treat the case of axially symmetric deformation taking the 
z-axis as the symmetry axis; /J^,- = ^ Pz- In this case, we can replace / dQ by 
jdcosO, Dijj^{Q) by di^^iO), R{Q) by Ry{e), and R{Q) by Ry{e). Here Ry{e) and 
Ry{9) are the rotation operator and matrix, respectively, representing the rotation 
by the angle around the y-axis. 

2.2. Elimination of spurious center- of -mass motion 

As already mentioned, the total center-of-mass motion is separated out of the 
wave function ^na{Px, Py- l^z)- However, the separation of the total center-of-mass 
motion is no more true for the angular-momentum-projected wave function 
^naif^x: Py, l^z)- Nevertheless the spurious admixture of the the total center-of-mass 
motion can be removed rather easily, as we explain now. 

By using Eq.(2.5) we separate the center-of-mass wave function from ^na{Px, Py, Pz) 



as 



- -^i^Gx + ^Gy) - -^^Gz] ^naiPx = Py,Pz), 



^ "2 

$na{l3x= Py,/3z) =A[ew{-Yl E -^i^ik - ^Ckf} (piai) 

1=1 k=x,y,z k 



(2-10) 



(2-11) 



where Q is some constant number. Clearly ^na{Px = /3y,Pz) does not contain 

the center-of-mass coordinate Xq and is the internal wave function of (Pnaifix = 
f3y,Pz)- The correct a-condenscd wave function with good angular momentum is 
not ^iaiPx = Py,(3z) given in Eq.(2.9) but is ^;(„(/3x = Py,l3z) defined as 

^iaWx = Py,Pz) = J dcOSedU{0)Ry{0) $naiPx = Py,Pz)- (2-12) 

The matrix element of a general translationally invariant scalar operator O with 
the correct CK-condensed wave function ^naWx = l3y,Pz) is calculated as 

{3UPx = l3y,Pz)\d\$L{Px = Py,l3z)) 
{$UPx=Py,Pz)\$LWx = Py,Pz)) 
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JdcOSedioie){$na{Px = Py,Pz)\Ry{9)\$na{Px = Py,Pz)) ' 



(2-13) 



The important technique for the calculation of the matrix elements of the integrands 
in this expression is to use the following relation 

{$na{Px = Py,Pz)\dRy{e)\$naiPx = Py,Pz)) 

= Q2^i'^na{Px = Py,Pz)\dRym^na{Px = Py,Pz)), (2-14) 

P(e) = (exp(- Y: §^^.)l^.W|exp(- E 

k=x,y,z ^ k=x,y,z ^ 



X 



\AnJ BlBl + {l/A){Bl-Blfs\n^e' 



(2-15) 



This relation holds because the operator O does not depend on the center-of-mass 
coordinate Xq- By using this relation we get 

{$-UP.=Py.Pz)\d\$L{Px=Py,Pz)) 
{$UP,=Py,P,)\$UPx = Py,Pz)) 

_ JdcOSed^o{9){^naiPx = Py,Pz)\ORy{0)Wna{Px = Py, [3,)) / P{9) 
J d COS edio{e){^na{Px = Py, Pz)\Ry{0)\^na{Px = Py,Pz))/P{0) ' 

(2-16) 

The explicit formulas of the matrix elements of the overlap, the kinetic energy, 
the two-body nuclear force, and the Coulomb force in the case of the ^Be system are 
given in the Appendix A. In Appendix B we will show that the integration over the 
Euler angle 9 can be performed analytically. 

2.3. Non-zero spin states on the Px{= Py) = Pz line 

On the Px{= py) = Pz line the intrinsic state ^na{Px = Py,Pz) is spherical 
and it may seem that we do not get non-zero spin states by the angular momentum 
projection. But we can construct the non-zero spin states even on the Px{= Py) = Pz 
line using angular momentum projection by the following usual limiting procedure: 
We first construct the non-zero spin state ^naiPx = PyiPz) of Eq.(2.12) for Px{= 
Py) Pz by angular momentum projection. Next we normalize this state obtaining 
the normalized state $^^{Px = Py,Pz), 

$^niiPx = Py,Pz) = ,|r!t^^!^iil ' (^-l^) 
\\^LiPx= Py,Pz)\\ 

where ||^|| denotes the norm of ||#|| = Since \\$naiPx = Py,Pz)\\ = 1) 

there exists the limiting state 

hm $^^{Px = Py,P,). (2-18) 
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This limiting state is the non-zero spin wave function we adopt on the Pxi= Py) = l^z 
hne. 

Since hm^g^^/j^ = I3y,(3z) = 0, the use of ^naif^x = fiy^Pz) near the 

hmit of (3x = (3z is not easy with respect to numerical accuracy. Therefore in our 
present study of the ^Be system, we calculated the matrix elements with ^2a{Px = 
f3y,Pz) analytically. This was made possible by performing the angular momentum 
projection analytically. Then the matrix elements with ^2a{Px = Py,Pz) near and 
on the l3x{= Py) = Pz line were obtained without any difficulty in the numerical 
calculations. Details are explained in Appendix C. 

2.4. Microscopic Hamiltonian 

The Hamiltonian H we use in this paper is written as 

H = T-Tg + Vn + Vc, (2.19) 

where T is the total kinetic energy operator, Tq is the center-of-mass kinetic enegy 
operator, Vn is the effective two-body nuclear force, and Vc is the Coulomb force 
between protons. Explicit forms of T and Tq are as follows 



1=1 ' ^ 



As the effective two-body nuclear force, we adopt the force No.l of VolkovS which 
has the following form, 

8 2 2 

^iv = ^E{(1 - M) - MP^P,},, ^ ^„exp(-^). (2-21) 

i^j n=X 

The Coulomb force is given as 

In our previous study of the alpha condensation in ^^C and ^^O, we used an effective 
nuclear force named Fl which contains a three-nucleon force in addition to the two- 
nucleon force eP. The reason for this was that in both nuclei we treated not only the 
a-condensed state of dilute density but also the normal-density states including the 
ground state. Therefore we had to use a nuclear force which has reasonable density- 
dependent properties. On the other hand, in the present study we treat only the 
ground band states of ^Be, all of which are expected to have a similar value of low 
density. Therefore we here use a simpler effective nuclear force than in the previous 
study. 



§3. Results 
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3.1. Energy surface of the intrinsic state 

We first calculate the binding energy for the spherical state of ^Be which has two 
parameters b and f3x = Py = Pz = Ro- In Fig-1 we give the contour map of the energy 
surface in the two parameter space, b and Px = (3y = (3z- The binding energy Ea 
of the single alpha particle which is independent of the Majorana parameter of the 
force No.l of Volkov takes its lowest value E^ = -27.08 MeV for the size parameter b 
= 1.37 fm. As the strength of the Majorana parameter M of the Volkov force in ^Be, 
we adopt M = 0.56. It is because with this value of M, the lowest binding energy 
takes the value -54.33 MeV as seen in Fig.l which is only about 0.17 MeV lower than 
the theoretical two-alpha threshold energy, 2Ea = —54.16 MeV, in agreement with 
the experimental binding energy from the two-alpha threshold energy, about + 0.1 
MeV. The value of the b parameter at the energy minimum is 6 = 1.36 fm which is 
very close to but slightly smallero than the b value of the free alpha particle. It is 
known that this combination of the values of b and M which are close to b= 1.37 
fm and M=0.56, respectively, gives a very good reproduction of the a-a scattering 
phase shifts within the framework of the resonating group methodcP . 

The qualitative feature of this map is similar to the corresponding maps in ^^C 
and ^^O given in our previous paperEl^. This map shows a valley running from the 
outer region with large value of f3x = Py = Pz > 10 fm and b ~ 1.36 fm. The valley 
has a saddle point at Px = Py = Pz ^ 8.5 fm. Beyond the saddle point the binding 
energy is almost equal to 2Ea. The height of the saddle point measured from the 
theoretical threshold energy, 2Ea, is about 0.56 MeV. The reason of the appearance 
of the saddle point is the same in the cases of ^^C and ^^O and is attributed to 
the increase of the Coulomb energy and kinetic energy towards the inward direction 
which around the saddle point region is not yet compensated by the gain in nuclear 
potential energy. The minimum of the energy surface is located at Px = Py = Pz ^ 
3 fm. 

The result that the a-condensed wave function reproduces well the experimental 
binding energy implies that the ground state of ®Be can be considered as a 2a-particle 
s-wave quasi-bound state with a " gas- like" structure of dilute density. By " gas-like" 
we mean that the state is weakly bound and spreads widely forming a dilute density 
state and constituent clusters move rather freely. 

By fixing the b parameter to b = 1.36 fm, which is the b value at the energy 
minimum of Fig.l, we now calculate the energy of the deformed states. Fig. 2 is the 
contour map of the energy surface in the two parameter space, Px{= Py) and Pz- This 
figure shows that the energy minimum lies on the Px = Pz hne, which means that 
the energy minimum is attained for the spherical shape. However, it is well known 
that nuclear deformation can not be deduced from a calculation without angular 
momentum projection. That is, we often encounter the case where the energy of 
the angular-momentum-projected state has a minimum for non-zero deformation 
although the energy of the intrinsic state without angular momentum projection 



* This smaller 6- value may be interpreted as an effect of the Pauli-principle:one way to minimise 
the repulsive effect of antisymmetrisation is that the clusters make theniselves smaller. This is clearly 
seen in the case of a low density gas of deuterons as explained in Ref-U 
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favors a spherical shape. 

3.2. Energy surfaces of spin-projected states 

In Fig. 3 we show the contour map of the energy surface corresponding to the 
spin-projected states J'^ = 0+ in the two parameter space, /3x{= Py) and P^. Now 
we see that the energy minimum no more hes on the l3x = (3z line but is located at 
Px{= Py) ~ 1-8 fm and ~ 7.8 fm in the prolate region of the map. The value of the 
energy minimum of -54.45 MeV means that the energy gain from angular momentum 
projection is only about 0.12 MeV. Reflecting this small energy gain, there is a valley 
with an almost fiat bottom running through this energy minimum to the second 
energy minimum at Pxi= Py) ~ 4.1 fm and Pz ~ 0.0 fm in the oblate region of the 
map. The bottom line of this valley crosses the Px = Pz line sX Px = Pz = fm 
which is just the energy minimum of Fig.2 for the spherical a-condensed state. The 
variation of the binding energy along the bottom line of the valley is less than 0.12 
MeV. 

How much different are the spin-projected wave functions along the bottom line 
of the valley from one another ? The squared overlap | l^n )P of the normalized 
wave function <Pj^ at the energy minimum and the normalized wave function at 
Px = Pz = ^ is about 0.98. This implies that is almost equivalent to ^n^. 
We also calculated the squared overlap of at the energy minimum 

and the normalized wave function ^^ii the second energy minimum. The value 
of Kf^j^ltpjj'f)!^ is about 0.99. This result implies that <Pj^ in the prolate region 
is almost equivalent to in the oblate region. Hence we can say that the 0+ 
component contained in the deformed intrinsic state along the bottom line of the 
valley is almost unchanged and is almost equivalent to the spherical condensed state 
&i Px = Pz = ^ fm which is the energy minimum of Fig.2. 

In Fig.4 (a) we show the squared overlap 0{px,Pz) = \{$2a'^ {Px, Pz)\$f)\'^ in 
the form of the contour map in the two parameter space, Pxi= Py) and P^, where 
$2a^{Px,Pz) is the normalized wave function of ^^iPx = Py,Pz) of Eq.(2.12). The 
feature of the contour map of Fig.4 (a) is rather similar to the contour map of Fig.3. 

Now we show in Fig. 5 the contour map of the energy surface for the spin- 
projected states to = 2"*" in the two parameter space, Px{= Py) and Pz- Unlike 
the 0+ case, we see no local minimum. The binding energy becomes deeper as Px or 
Pz becomes larger. Let us consider the region between two lines, /J^; + « 3 f m and 
Px + Pz ^ fni. The line, Px + Pz ~ ^ fm, is just the beginning of the energetically 
steeply-rising region towards Px = /S^ = 0. The line, Px + Pz ~ 7 fm, is approximately 
the boarder of the interaction region of two alpha particles. The region between the 
two lines, /3a; -|- ~ 3 f m and Px + Pz ^ is rather flat in energy, that is the bind- 
ing energy in this region lies approximately between -51.0 MeV and -51.8 MeV. The 
height of this plateau region measured from the 0^ energy minimum, -54.45 MeV, 
is approximately between 2.65 MeV and 3.45 MeV. Although, without imposing the 
correct boundary condition of the resonance, we cannot make a definite statement, 
we can safely conjecture that we will have a 2^ resonance of the excitation energy 
between 2.65 MeV and 3.45 MeV. This excitation energy is in good agreement with 
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the experimental 2+ excitation energy at 2.9 MeV. The reason of our conjecture is 
as follows. Let us make the strength of the Majorana exchange mixture M slightly 
smaller, taking M= 0.56 instead of M= 0.54. The contour map of the J'^ = 2"*" en- 
ergy surface for M = 0.54 is shown in Fig. 6. The feature of the contour map of Fig. 6 
is very similar to that of Fig. 5, but we now have a local minimum at Px{= Py) ^ 
3.9 fm and /J^ 0.0 fm in the oblate region of the map. The binding energy of this 
local minimum is -52.28 MeV. We show in Fig. 7 the contour map of the J'^ = 0^ 
energy surface for the same M = 0.54. The feature of the contour map of Fig. 7 
is very similar to that of Fig. 3 with M = 0.56 and the energy minimum which is 
slightly deeper than that of Fig.3 is now at -55.19 MeV. Then the excitation energy 
of the 2+ state is 2.91 MeV. This value of the 2+ excitation energy is very close to 
the experimental excitation energy of 2.9 MeV. 

In Fig.4 (b) we show the squared overlap 0{Px,Pz) = \{$2a^ iPx, Pz)\$i^ )\'^ in 
the form of the contour map in the two parameter space, Px{= Py) and Pz, where 
02a^ {Pxj Pz) is the normalized wave function of $l'^{Px = Py,pz) of Eq.(2.12), and 

is the normalized wave function of the energy minimum of Fig. 6, Px{= Py) ^ 
3.9 fm and Pz ~ 0.0 fm. This figure shows that the 2+ wave function is almost 
unchanged from along the valley running from this energy minimum in the 
oblate region deeply into the region of prolate deformation. Hence we can say that 
the 2"*' component contained in the deformed intrinsic state along this bottom line 
of the valley is almost unchanged. 

We also checked the cases of M < 0.54 for the 2+-state and again found the 
existence of a local minimum. The excitation energy of the 2^ state at the local 
minimum was found to be around 2.9 MeV. 

From the above results for the 0+ and 2+ states, in using the deformed 2a-state, 
we reach two conclusions: (1) For the 0"*" state, the introduction of the deformed 
2a-state plays essentially no role and the description with the spherical 2a-state 
is justified. (2) The deformed 2a;-state yields a nice description of the 2+ state, 
which gives us a new understanding about the character of the 2+ state. Namely 
the 2+ state can be considered to be a gas-like state of 2Q:-particles in a relative 
d-wave. The merit of the introduction of the deformed 2a-state is just this point. 
That is, we have obtained the picture about ^Be that not only the ground state 
but also the excited 2+ state can be interpreted in terms of a gas-like 2a-structure. 
This conclusion encourages us to introduce the deformed wave function of the a- 
condensation in ^^C and heavier self-conjugate 4n nuclei for a future study of the 
non-zero spin excitations of a-condensate states. 

Let us now study the 4+ state. Fig.8 shows the contour map of the energy surface 
corresponding to the spin-projected states of J'^ = 4+ in the two parameter space, 
Px{= Py) and [3z- Like the 2^ case there is no local minimum. We checked the cases of 
smaller M values than 0.56 including 0.54 but we found no local minimun in all cases. 
In Fig.8 we see that the approximately flat region of the 4+ surface is slightly pulled 
into a smaller f3x + Pz region compared to the case of the 2^ energy surface. The 
binding energy in the region between the two lines, Px+Pz ~ 2 fm and Px+Pz ~ 4 fm, 
lies approximately between -42 MeV and -45 MeV, which in terms of the excitation 
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energy is between 9.45 MeV and 12.45 MeV. Although this excitation energy region 
between 9.45 MeV and 12.45 MeV is near the experimantal 4+ excitation energy, 11.4 
MeV, wc need more careful calculation including the resonance boundary condition 
in order to get reliable conclusion about the 4"'" resonance. 

In this context it is useful to mention the calculation of the energy curve of 
^Be using Brink's 2-a wave function. The energy curve is the binding energy as a 
function of the inter-ct distance R. The energy curve for the 2"*" resonance has a 
local minimum inside the interaction region with R < 6 fm, while that for the 4+ 
resonance does not have any local minimum. 

3.3. Comparison with the solution of the 2a Hill-Wheeler equation 

The wave function of the 2a-condensed state belongs to the functional space 
spanned by the wave functions of the following form, 

A{xir)<P{a,)(l>{a2)}, (3-1) 

with x(r) representing an arbitrary function. Here r is the relative coordinate be- 
tween two a clusters, r = Xi — X2. The diagonalization of the Hamiltonian within 
this functional space is accomplished just by solving the well-known equation of the 
resonating group method (RGM), 

{<P{a,)(P{a2)\{H - E)\A{x{r)<P{a,)<p{a2)}) = 0. (3-2) 

The solution of this RGM equation is equivalent to the solution of the 2a Hill- Wheeler 
equation, 

J d'R -ls)\{H - E)\^^{^R, -Ifi)) f{R) = 0. (3-3) 

Our wave functions of 2a-condensation are approximations to the RGM wave func- 
tions A{x{f)4'{'^i)'P{(^2)} obtained by solving Eq.(3.2) in imposing wave functions of 
Gaussian shapes. Since the numerical solution of the 2a RGM equation of Eq.(3.2) 
is not difficult with present-day computers, one may ask why we discuss approxima- 
tions to the RGM wave functions. The answer is of course obvious: Our study in 
this paper is to extract and elucidate the a-condensation character of the '^Be states 
which the RGM study of ^Be until now has not clarified even though it has given 
numerical values of the ^Be wave functions. Furthermore our present study of the 
^Be states is made in the scope which is not restricted only to the ^Be states but ex- 
tends to excited states of general self-conjugate 4n nuclei which are of dilute density 
and are expected to exist near the na breakup threshold. For the self-conjugate 4n 
nuclei, heavier than ^Be, especially for their excited states with dilute density, it is 
no longer easy but very difficult, especially for n > 4, to solve the RGM equation of 
motion, 

(0(ai) • • • cl>ian)\iH - ^)|^{x(^i, • • • €„-i)0(ai) • • • = 0, (3-4) 

where ^1, • • • are the relative coordinates between na clusters. 
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For 0"'", we solved the 2a Hill- Wheeler equation of Eq.(3.3) and compared the 
solution with our wave function The Hill- Wheeler equation was solved by 

discretizing the radial integration as follows 

- E)\^^'^{Rj)) f{Rj) = 0, 

j 

^^'\r)^ j<fR^'^{\R,-\R), (3-5) 

where R is the polar angle of R. The energy expectation value of as 
a function of R, {<P^^^ {R)\H\^^^^ (R)) / {^^^^ {R)\<P^^^ (R)), takes its minimum at 
R = 3.45 fm with the minimum energy -52.069 MeV. It should be noted that this 
value is higher than the energy (-54.448 MeV) of by 2.379 MeV, which is a rather 
large value. By adopting Rj = 0.5 x j fm with j = 1 ~ 23, we obtained -54.444 
MeV as the lowest eigen energy, and the squared overlap of the corresponding wave 
function with our wave function (pj^ is 0.9973. When we adopt one more mesh point 
as Rj = 0.5 X j fm with j = 1 ~ 24, we obtained -54.446 MeV as the lowest eigen 
energy, and the squared overlap of the corresponding wave function with our wave 
function <Pj^ becomes 0.9980. We can say that we have obtained almost converged 
value for the lowest eigen energy of the Hill- Wheeler equation. In principle, the 
lowest eigen energy of the Hill- Wheeler equation of Eq.(3.5) should be lower than 
the energy of . But the above values of the lowest eigen energy are still slightly 
higher than the energy of . Prom this result we can say that our wave function 
is almost exactly equal to the wave function of the lowest energy solution of the 
Hill- Wheeler equation. This conclusion is surprising but it clearly shows that our 
model wave function ^2a^ iPx = f^yif^z) is very much suited to the ground state of 
^Be which is a threshold state with dilute density. We can say that this conclusion 
strongly supports our new picture of the ^Be structure as a very dilute gas-like 
structure of 2 a-particles rather than a dumb-bell structure of 2 a-particles. 

§4. Discussion 

Our study in this paper is based on the use of the spin-projected state ^2a{Px = 
Py,Pz) obtained from the deformed intrinsic state ^2a{Px = Py,Pz)- Our deformed 
intrinsic state ^2a{Px = Py, Pz) ^s, characterized as having a gas-like 2Q!-structure,i.e. 
a weakly bound 2a-state. However, in order to be able to impose the same character 
to the spin-projected state ^2aWx = f^y^Pz)-, it is necessary that the spin-projected 
state ^^ail^x = l3y,l3z) is contained in the intrinsic state 'l'2a{(^x = l3y,(3z) with 
non-negligible amplitude, where ^^^if^x = Py^l^z) and ^2aiPx = Py,Pz) are the 
normalized states of ^iaif^x = Py,Pz) and ^2a{Px = Py,l3z), respectively. Therefore 
we have calculated the magnitude of the squared overlap amplitude, \{^2ail^x = 
Py,(3z)\$2a{Px = l^y, Pz))? , for J = ~ 4, and the results are given in Figs. 9 (a) ~ 
(c) in the form of the contour map in the two parameter space, (ix{= f^y) and (3z- 
We see that the minimum-energy states and have strong overlap with their 
respective instrinsic states. 
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The energy surfaces given in Figs. 3 and 5 are fairly flat. But at the same time 
we see in Figs.4 (a) and (b) that the wave functions in the energetically flat region 
are very similar to one another. Therefore the flatness of the energy surface does 
not imply the appearance of the low-energy excited states. Also the superposition of 
the wave functions by solving the Hill- Wheeler equation will result in only a slight 
change of the energy and wave function from the energy and wave function at the 
energy minimum of the energy surface, respectively. 

As explained in § 2.3, the non-zero spin states on the Px = Pz line were con- 
structed by the limiting procedure and the matrix elements of these limiting states 
were calculated analytically. As is seen in the maps of Figs. 5, 6, and 8, the contour 
lines cross the (3x = Pz line smoothly with no singular behavior near the Px = Pz 
line. In Fig. 10 we show the binding energy of the 2+ state along the the Px = Pz 
line in the case of M = 0.56. 

In Fig.4(a) we saw that the minimum-energy 0+ wave function which is projected 
out from the prolate intrinsic state is almost the same as the second- minimum-energy 
0"'" wave function which is projected out from the oblate intrinsic state. Similarly 
we saw in Fig.4(b) that the 2+ wave function which gives the local energy minimum 
and which is projected out from the oblate intrinsic state is almost the same as 
the 2"*" wave functions which are projected out from the prolate intrinsic states. 
These results look very strange at first sight. In order to check whether the prolate 
intrinsic state is very different from the oblate intrinsic state or not, we calculated 
the overlap between the intrinsic wave function at the energy minimum for 0"*" 
which has prolate deformation with the general intrinsic state ^2a{Px = Py,Pz) 
including the intrinsic wave function of the second energy minimum which has oblate 
deformation. In Fig. 11 we show the squared overlap \{i>2a{Px, Pz)\^f)\'^ in the form 
of a contour map in the two parameter space, Pxi= Py) and Pz-, where ^^^{Px,Pz) 
is the normalized wave function of ^2a{Px = Py, Pz) and is the normalized wave 
function of (pj. We see that the overlaps between the prolate intrinsic state and 
the oblate intrinsic states near the second energy minimum are very small. Below, we 
explain the reason why we have such seemingly strange results for the spin-projected 
states and indicate that this is due to the simplicity or the smallness of the functional 
space of the 2a system. The intrinsic wave function of our model which has the form 
of Eq.(2.11) with n = 2 is written as 

$2a{Px = Py,Pz)=A{exp(^-^-^!^^ - ^) </,(ai)</,(a2) } , 

rk = Xik-X2k, {k = x,y,z). (4-1) 



The relative wave function exp{— (r^ -l-r^)/i?^ — r^/B^} can be expanded as follows 




+ •••, (4-2) 
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where f is the polar angle of r = Xi — X2- The sign of the coefficient of the 
second term, (l/SSf — 1/35^), is opposite for prolate and oblate deformations, and 
it makes the overlap between prolate and oblate intrinsic states small. However once 
we project out the good spin state, the situation changes strongly. For the 0"'"-state, 
the projected wave function has the same following form for the leading term 



irrespectively of the deformation of its original intrinsic state. Similarly for 2"*", the 
projected wave function has the same following form for the leading term 



irrespectively of the deformation of its original intrinsic state. Therefore the overlap 
between spin-projected wave functions cannot become small but rather becomes close 
to unity as long as the values of {2/B'^ + I/BI) are close to one another, irrespectively 
of the deformations of their intrinsic states. We can say that the essence of our above 
analysis is as follows: In the 2a system the spin-projected wave function is governed 
by a Gaussian-packet-like function Xj{f) of a single coordinate r, the relative distance 
of the two a's. Thus as long as one global character such as (xj^^Ixj) is similar, 
the spin-projected wave function can not differ so much, and there is no room for a 
difference with the original deformation. We see in summary that the origin of the 
seemingly strange results about the spin-projected wave functions comes from the 
simplicity or the smallness of the functional space of the 2a system. Hence this kind 
of things will not be possible to occur for 3q; and heavier systems. 



In order to study non-zero-spin excitations of a-cluster condensation in the self 
conjugate 4n nuclei which we proposed in our previous paper, we introduced a spa- 
tial deformation into the a-cluster condensate. The wave function of the a-cluster 
condensate with deformation is obtained by a slight and natural modification of our 
previous wave function of the spherical condensate of a-clusters. As a first step and 
test case for the study of deformation, we, in this paper, investigated the ^Be nucleus 
in applying the wave function of the deformed a-cluster condensate to the 2a-system. 
The 0+ ground state, the 2+ excited state at 2.9 MeV with = 1.45 Mev, and the 
4+ excited state at 11.4 MeV with ~ 7 MeV are known to constitute a rotational 
band whose intrinsic state has a 2a cluster structure with large inter-a distance. Our 
study in this paper gave us the following conclusions: (1) the spherical a-condensed 
wave function reproduces well the binding energy of the ground state, which implies 
that the 0+ ground state can be considered as a 2a-particle quasi-bound state in 
relative s-wave with gas-like structure of dilute density. (2) The extension of the 
model space so as to include the deformed structure gives us only a slight energy 
gain of about 0.12 MeV for the 0"*" ground state. Corresponding to this small change 
of the binding energy, the change of the wave function due to the extension of the 




(4-3) 




(4-4) 
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model space is very small, although the deformed intrinsic wave function before the 
spin projection to J = at the energy minimum is quite different from the spherical 
optimum wave function. Thus as far as the O"*" ground state is concerned, the intro- 
duction of the deformed state plays essentially no role. (3) Our 0+ wave function is 
found to be almost exactly equal to the 0+ wave function obtained by the generator 
coordinate method using Brink's 2a wave function. This result strongly supports 
our new picture of the ^Be structure that it is more a very dilute gas-like structure 
of 2Q;-particles rather than a dumb-bell structure of 2a-particles. (4) The deformed 
wave function reproduces well the observed excitation energy at 2.9 MeV of the 2+ 
state. This result gives us the picture that not only the ground state but also the 
excited 2+ state of ^Be can be considered as a 2a-state with a gas-like structure 
of dilute density. This conclusion encourages us to introduce the deformed wave 
function of the a-condensation also for ^■^C and heavier self-conjugate 4n nuclei for 
a future study of non-zero spin excitations of the a-condensation. 

It is worth mentioning that our a-particlc condensed wave function, introduced 
in Ref.[l], has so far reproduced quite accurately 0+ -threshold states in ^^C and ^^O 
and the rotational states in ^Be, treated in this paper. The fact that these results are 
obtained without a single adjustable parameter is, we believe, quite remarkable and 
gives strong credit to the correctness of the physics contained in our wave function. 
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Appendix A 

Matrix elements 



The deformed condensed wave function of Eq.(2.6) is written for the 2a system 



as 



^2a{(3x=Py,l3z) = J Ri d^R2 exp(-^ Yl '^''iRi,R2). (A-1) 

1=1 k=x,y,z k 

By transforming Ri and R2 into the center-of-mass vector Rq and the relative 
vector ii, 

Ri = Rg + \r, R2 = Rg-\r, (A-2) 
this wave function is rewritten as 

^2a(/?. = Py.Pz) = 4! (^)'^' / d'Ro exp { - J2 llRlk - ^{Xg - Rg?) 

k=x^y,z ' ^ 

/I 1 
d^E exp(- ^^fc) -4[exp{-^(r-i«)2}<^(ai)0(a2)], 

k=x,y,z ' k 

where Xq = {Xi + ^2)/2 and r = Xi — X2, with Xi denoting the center-of- 
mass coordinate of the i-th a cluster, a^. Therefore the internal wave function, 
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$2aii3x = Py,/3z), of Eq.(2.11) is written as 



$2a{Px = Py,Pz) = Q' J d^R exp(- ^1 J^^k) 



k=x,y,z 



xA 



exp I - ^(r - Rf^<p{ai)<^{a2) 



Now we introduce the following wave function, ^2a{Px = l3y,Pz)-, 

^2a{Px = Py,Pz) = 4! (^)'^'^ exp(-lx^)^2a(/9:. = Py,Pz 



K=x,y,z « 

It is easy to prove the following relation 

i^LiPx = ^y,Pz)\d\$L{Px = Py,f3z)) 
{HaWx = Py,Pz)\$LiPx = Py,Pz)) 

^ Jdcose4o{e){^2a{(3x = Py,Pz)\ORy{0)W2a{Px = Py,Pz)) 
jdcOSedUe){W2a{f3x = f3y,Pz)\Ry{0W2a{f3x = Py,Pz)) ' 

Below we give the explicit formulas of the matrix elements of the overlap, the 
kinetic energy, the two-body nuclear force, and the Coulomb force in the 2a system. 
We use below the notation 



(A-3) 



The formula of the overlap is as follows. 

{^2a{l3x = Py,l3z)\Ry{9)\^2aiPx = l3y,Pz)) 



= I d'Sd'R exp {- E Sl + {Ryi9)R)l ^ ^^B(^) j ^B(^^^ 

k=x,y,z '~^k 

= i27rfPtf3! (G0 + G1 + G2), 

2 



Gn = 



Gi = - 



(2a? + l){(2a! + l)(2ai + l)+7} 



^(iaf + l)(iaf + l){(iaf + !)(§«? + l)(iai + l)(|ai + 1) + i} ' 



Go 



(al + inal + l) ■ 
In the above relation, we introduced variables, ai, a2, and 7, defined as 

ai = (3^/b = py/b, a2 = f3z/b, 7 = (a?-ai)2sin2 0. (A-4) 
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The formula of the kinetic energy is as fohows. 

{^2aiPx = Py,f3,)\{T-TG)Ry{e)\^2a{Px = Py,Pz)) 



J ^ 1. ^ 2/5/, J 



k=x,y,z 



4 



Saj + l al + al + 1 



2a? + 1 {2al + l){(2a2 + l){2al + 1) + 7} 



{Sajal - 4)(f afai + af + + 1) 



4(|af + 3al + 1) 

(ia? + l)(|a? + l)(lai + l)(iai + l) + | ' (^af + l)(|af + 1) 



+ 



G2 2af(ai + l)+ai(af + l) \ 



The formula of the two-body nuclear force is as follows. 

{^2a{Px = Py,Pz)\VNRym^2a{Px = Py,Pz)) 

= / d'Sd'Rexp{- Y: ^^±^^^} {$^{S)\Vr,\$^{R)) 

k=x,y,z '"^^ 
n=l " 



4Xd 



{2al + 1) Ka? + l){(2a2 + 1) + 1) (2a2 + 1) (^.^^i + 1) + K - 2)^ ^1 

4X, 



(ia2 + l){K + |)a2 + l}[(lc,2 + l)|(^^+3)^2 + i|(1^2 + i)|(^^ + 3)^2 + i| + (^^ + l^ 

4(Xe - 2Xd) 



(§af + l){K + i)a? + l}[(|a? + l){K + i)a? + l}(iai + l){K + i)ai + l}+K-l)2| 

16(Xrf + Xe) 



(a2+l){(in„+l)af+l}-^l [(a2+i)|( i^^+i)c,2+i|_^^l [(«2+i)|^i„^+i^^2+i|_^ 
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4{Xd - 2Xe) 



(a2 + l){K + l)a2 + l}[(a2 + l){K + l)a2 + i}(a2 + i){(^„ + l)a2 + i|+^2 2 
16{Xd + X,) 



(af + l){(iu„+l)af + l}^(ai + l){(iu„+l)ai + l}/ ' 

As explained in section §2.4, since we make use of the Volkov interaction as the 
two-body force, we should note that Vn reads as follows, 

8 2 2 

i^i n=l '"^ 

In the above formula we used the variables defined below, 

= 8 - lOM, Xe = lOM - 2, 
62 



a2 + 262 ■ 



As for the Coulomb force, we express it as a superposition of Gaussian functions 
which takes the form, 

1^3 V 



2^^ —V^Jo ^^P^-^^^-^ )■ 



Therfore, we can easily obtain the matrix elements by changing the variables used 
in Vjv in the following way, 

2 1 

^/"Va2 + 262J 0Fio Vi + 262^2j 

It is convenient to change the integration variable from 77 to a new variable ^ defined 
as 

^2 262j72 n~ j 1 s3/2 

The range of integration over ^ is from zero to one. The numerical integration over 
^ was performed by using the Gauss-Legendre method. 
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Appendix B 

Analytical calculation of spin projection 

Matrix elements with the angular- momentum-projected wave functions are ob- 
tained by integrating the formulas of Appendix A multiplied by dQQ{6) over the angle 
9. We should note that in the formulas of Appendix A only 7 = (af — a^)^ sin^ 9 
depends on the angle 9. When limited to the case of 0+, 2+, and 4+, d^o{9) contains 
only cos^ 9 and cos*^ 9 as the terms which depend on 9. Therefore, the integrals for 
0"*", 2+, and 4+ are classified into the following three types, 



/I -j^ 1 
dx , = 2arcsin(— -), 

VA'^ - 1 + A'^ arcsin(^^) , 



1 ^2 
dx 



1 VA^ - x2 



1 3,4 
dx- 



4(2 + 3^2) - 1 + arcsin(^) , (B-1) 



in the case of the overlap and nuclear force terms, and into the following three types. 



dx- 



1 (^2 _ 3,2)3/2 ^2^/]423l' 



^ dx—— — ■ — = , — 2arcsinf— V 

-1 (A2-x2)3/2 y'^2_ 1 V^^' 
/■I / 1 \ 

/ dx— -— = V^2_i^ =-3^^arcsin - , (B-2) 

y_i (^2 _ 3.2)3/2 _ I \Ar 

in the case of the kinetic terms. Here x = cos 9 and we note that A always takes 
values more than unity. 

Appendix C 

Limit of Px Pz for spin projected states 



We stated in §2.3, that there exist non-zero spin states on the Px{= Py) = Pz line. 
They are defined by the limiting procedure of Px Pz- We can obtain the matrix 
elements with these limiting states analytically by performing Taylor expansion with 
respect to s = af _ q,2 around s = 0. In this appendix, we display the series 
expansions with respect to s up to the second order for the case of the 2"*" state. For 
the sake of the series expansion, we introduce i = -|- and rewrite and 02 by 
s and t as 

al = lis + t), al = ^it-s). (C-1) 



By using 



^oo(^) = - ^, (C-2) 



we get the series expansions about the overlap terms as follows, 

/-i /3 2 1\ ^ 4 1 2 2 1 3 
Go = J_^dx(-x -^)Go = j-^^^j^s -j-5^^jY)^s + 
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15{(it + l)(|t + l)}V2' 



30{(it+l)(ft + l)}9/2' 



Go 



dx{ 



^2 2' 
About the kinetic terms, we get 



1 



- ^ {^2a{Px = Py, PzW - TG)Ry{x)\^2aiPx = Py, Pz)) 



(C-3) 
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(Go + Gi + G2) 



1 3t-4 



1 2 



15(^ + 1)9/2 10 (t + 1)11/2 

1 15*2 + 24t - 32 



, 1 15t''^ + 56*2 _ ;^28 „ 
rS - — ^T^^ — T -ttttttS-^ + • • • 



120 {(I* + l){\t + l)}9/2 ■ 640 {(1* + l){\t + 1)}1V2 • 



-0 



About the terms which come from the nuclear force, we have 

j^J^x" - ]^{^2a{Px = Py,Pz)\VNRy{x)W2a{Px = Py,Pz)) 



(C-4) 



< \3/2 



(2.)Wi;.„(;^ 



n=l ^ 



4(X. + X.)(|44) 

, 4 V2(n„-2)2 2 2 y2(u„-2)2(2^x„t + n + 2) 3 



15 



{{t+l){Unt + 2)] 



7/2" 15 



9/2 



+4Xe 



(i^n + 1)' 



30 



{(|i + l)(^t + i)} 



7/2 120 (n 



{{t + l)iUnt + 2)} 

1 (n„ + l)2(2«^i + ^,„ + 2 

{(i*+l)(Mt+l)}'/' 



+4(Xe - 2Xd) 



-16(Xd + Xe) 



(Un - If 



30 



{(f* + l)(^i + i)} 
1 1 



V2 120 
1 



1 (n„-l)2(Mt + ^^ + 2] 
t + l)(^H^* + l)}'^' 



■s' + 



30 



4 '^"•"2"'"^ ^3_^ 



(Mt2 + + 1) 120 |^2«^^2 + + 1) 



9/2' 
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64 



2 64 {{Un + l)t + Un + 2}ul 3 



15 



(t + 2){K + l)t + 2} 



7/2^ 15 r 



{t + 2){{Un + l)t + 2} 



9/2 



+16iXd + Xe) X ^ . 



(C-5) 
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Fig.l: Contour map of the energy surface of the spherical state for M — 0.56 in the 
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Fig. 4(a): Contour map of the squared overlap 0{f3x, Pz) between the 0^ wave func- 
tion with (3x = Py ^ 1-8 fm, (3z = 7.8 fm and the 0+ wave function with variable 
Px(.= /3y) and /3z- Numbers attached to the contour lines are squared overlap values. 
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Fig. 4(b): Contour map of the squared overlap 0{(3x, Pz) between the 2+ wave func- 
tion with /3x = Py = 3.9 fm, /J^ = 0.0 fm and the 2+ wave function with variable 
l3x{= 0y) and Numbers attached to the contour lines are squared overlap values. 
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Fig. 5: Contour map of the energy surface of the 2+ state for M = 0.56 in the two 
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Fig. 8: Contour map of the energy surface of the 4+ state for M = 0.56 in the two 
parameter space, Pxi= Py) and (3z- 
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Fig. 10: Binding energy of the 2+ state on the (3x = (3y = (3z line. 
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